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Abstract 

We  give  a  simple  review  of  some  of  the  theoretical  aspects  of  Stimulated  Raman  Scattering 
with  particular  emphasis  on  those  features  which  would  be  of  experimental  interest. 

Ever  since  Chu  and  Scott^  first  demonstrated  that  the  Stimulated  Raman  Scattering  (SRS) 
equations  were  integrable,  these  equations  have  proven  themselves  to  be  quite  interesting  from  a 
theoretical  point  of  view'.  These  equations  have  been  intensely  studied^"®  theoretically  and  muneri- 
cally.  We  shall  discuss  here  some  of  th'  more  elementary  and  interesting  aspects  of  these  equations, 
with  particular  emphasis  on  those  simple  features  that  would  be  observed  experimentially.  W'e  shaiU 
avoid  as  much  as  possible,  the  complexities  of  the  Inverse  Scattering  Transform  (IST)  and  restrict 
our  discussions  to  those  points  that  experimentalists  would  be  interested  in. 

There  are  two  forms  of  these  equations.  The  first  is  the  more  genercd  case  where  it  is  possible 
for  the  ground  state  of  the  molecules  '■ould  become  depleted^'^.  This  case,  at  the  present,  cannot 
be  achieved  experimeutally,  so  we  shall  not  consider  or  discuss  it  here.  The  usual  experimential 
situation  is  where  only  a  very  smcdl  fraction  of  the  available  molecules  can  be  excited.  This  is 
the  case  originally  studied  by  Scott  zmd  Chu^  and  is  the  only  case  that  we  shall  consider  here. 
These  reduced  equations  are  what  one  usu2dly  czdls  the  SRS  equations.  There  have  been  some  very 
enlightening  numerical  studies®"®  of  these  equations  in  addition  to  all  the  analytical  studies.  And 
there  have  been  very  interesting  zmalytical  and  experimental  studies  of  SRS  soliton  evolution  in  the 
presence  of  damping®'®.  The  standard  form  of  the  reduced  SRS  equations  are 


=  —XA2 


(la) 
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d^A2  =  X*Ai 


(16) 


drX-hyX  =  AiA*2  (Ic) 

where  we  have  allowed  for  a  phenomenological  molecular  damping  factor,  7.  For  pulses  with  a 
sufficiently  short  width,  one  can  neglect  this  damping  factor.  We  shall  need  to  do  this  when  we 
discuss  the  soliton  solutions  but  for  now,  we  shall  retain  this  term.  This  damping  factor  will  be  used 
only  in  the  first  part  of  our  discussions,  after  which  we  shall  set  it  equal  to  zero.  In  (1),  Ai{A2)  is 
the  envelope  of  the  pump  (Stokes)  wave  and  X  is  the  material  excitation,  x  ^d  r  are  the  comoving 
coordinates 


X  =  z  (2a) 

T  =  t  —  z/u  (26) 

where  2  is  the  spatial  coordinate,  t  is  the  laboratory  time  and  u  is  the  common  group  velocity  of 
the  pump  and  Stokes  envelopes. 

The  nonlinearities  in  (1)  are  quadratic  nonlinearities  and  there  is  no  linear  dispersion.  Thus 
these  equations  are  identical  in  form  to  the  standard  3-wave  parametric  equations,  called  the  3- 
wave  resonant  interaction^®  (3WR1).  However  the  SRS  equations  are  a  degenerate  case  of  these 
equations  since  the  group  velocities  of  the  ptimp  and  the  Stokes  Eu-e  usually  taken  to  be  equail. 
Because  of  this,  the  SRS  equations  do  not  fit  into  the  stamdard  classification  scheme^®  of  the  3WRI 
equation,  but  rather  they  fedl  between  two  of  the  3WR1  cases.  The  classification  scheme  for  the 
3WRI  cases  depends  on  the  ordering  of  the  group  velocities  of  the  three  waves  and  whether  the 
pump  wave  (the  high  frequency  wave)  lies  between  the  other  two  group  velocities  or  not.  If  it  does 
lie  between  them,  then  we  have  the  the  soltion  decay  case  and  if  it  does  not,  then  we  have  the 
stimulated  backscatter  case.  With  the  above,  it  is  rather  easy  to  see  that  the  SRS  equations  are  a 
degenerate  case.  Just  note  that  if  the  group  velocity  of  the  Stokes  wave  were  slightly  higher  than 
that  of  the  pump,  then  the  pump  would  have  the  middle  group  velocity  and  the  SRS  equations 
could  then  be  classified  as  the  soliton  decay  case.  However  if  we  just  sfightly  shifted  the  group 
velocity  of  the  Stokes  pulse  to  be  just  below  that  of  the  pump,  then  we  would  have  the  other  case 
of  stimulated  backscatter.  Of  comse,  in  actuality,  Einy  experiment  would  fit  into  one  or  the  other  of 


2 


these  case.  But  also  in  any  experiment,  the  interaction  length  would  never  be  long  enough  for  the 
extremely  small  difference  in  the  group  velocities  to  have  any  effect.  Thus  in  effect  what  we  have 
is  really  a  hybrid  of  the  featiu’es  of  both  of  these  two  fundaunentcd  three-wave  interactions.  This 
demonstrates  itself  in  that  we  see  decay  of  the  pump  into  daughter  waves  (the  dominate  feature 
of  soliton  decay)  zmd  also  growth  of  nonlinear  oscillations  (the  dominate  feature  of  the  stimulated 
backscatter  case).  On  the  other  hand,  some  of  the  standcird  features  of  the  standard  cases  will  be 
absent.  For  excimple,  the  soliton  decay  case  always  has  a  threshhold  area  in  order  for  decay  to 
occur.  However  due  to  the  approximate  equality  of  the  group  velocities,  there  is  no  threshhold  area 
for  decay  of  the  piunp  in  SRS  because  the  envelopes  copropagate  and  never  separate.  (In  soliton 
decay,  the  threshhold  area  can  be  interpreted  as  the  width  necessary  in  order  for  the  interaction  to 
build  up  before  the  envelopes  separate.) 

One  of  the  main  features  of  the  SRS  equations  is  the  decay  of  the  pump  into  a  Stokes  wave. 
This  starts  as  a  linear  Raman  instability  and  can  best  be  seen  by  considering  (1)  in  the  hmit  of  very 
slowly  varjing  (in  time)  envelopes.  In  this  c<ise,  we  can  ignore  the  time  derivative  in  (Ic)  since  the 
damping  term  dominates.  Then  we  can  then  solve  directly  for  X,  obtaining  X  =  AiA^/j.  From 
(la)  and  (lb)  we  obtain 


(3a) 


- (^AyAi)A2  —  0 
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(36) 


Clearly  (3a)  is  decaying  while  (3b)  is  growing.  Consequently,  the  Stokes  wave  will  growm  from  any 
noise  at  the  expense  of  the  pump. 

In  general,  the  stable  configuration  will  always  be  a  Stokes  wave  with  no  pump  wave  remaining. 
It  is  interesting  to  note  also  that  it  is  only  on  such  a  backgroimd  that  one  can  have  a  soliton  solution. 
Of  course,  solitons  do  not  generally  exist  in  the  presence  of  any  damping,  so  let  us  now  assume 
that  the  time  derivative  will  be  so  large  (short  pulse  widths)  that  the  damping  constant,  7,  can  be 
neglected.  Then  a  soliton  solution  of  (1)  is 
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A{T)e' 
cosh  Z 
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(4a) 


By . . . 
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A2  =  A{t)  tanh  Z 


(46) 


(4c) 


X  =  ^— 
coshZ 

where  the  arguement  of  the  hyperbolic  functions  is 

Z  =  2itx-^f_^A\s)d$  (5) 

In  the  above,  A{t)  is  a  real  function  of  t  (otherwise  axbitrauy)  and  t)  and  <f>  are  real  constants.  This 
solution  represents  a  localized  disturbance  on  a  Stokes  bcickground  wherein  the  Stokes  ampUtude 
passes  through  zero  and  changes  sign.  This  solution  also  has  the  feature  that 


AlAi  +  AlA3  =  A^T)  (6) 

which  is  nothing  more  than  the  conservation  of  photon  number.  (From  the  derivation  of  (1),  one 
can  determine  that  the  squared  magnitude  of  the  seeded  amplitudes  An  are  indeed  photon  numbers 
fluxes.)  For  every  pump  photon  that  is  lost,  a  Stokes  photon  is  created  and  vice  versa.  One  can 
then  interpret  A^{t)  as  the  total  photon  number  flux  as  a  function  of  the  retarded  time.  If  we  fix 
X,  then  A‘^{t)  is  simply  the  total  flux  of  photons  seen  at  x  as  a  function  of  t.  And  as  the  pump 
and  Stokes  envelopes  copropagate  through  the  medium,  this  profile  of  A^{t)  never  changes  simply 
because  for  every  photon  lost  in  one  envelope,  another  one  has  appeared  in  the  other. 

Another  feature  of  these  soliton  solutions  is  that  they  are  in  general  transient  solutions.  Given 
any  pulse  of  finite  duration,  eventually  any  soliton  solution  will  vanish.  To  see  this,  let  us  take  an 
initial  pulse  profile  where  A(t)^  has  a  finite  area  and  also  has  some  knid  of  a  slowly  decaying  tail 
structure.  The  exact  structure  is  not  importemt,  only  that  it  eventually  vanishes.  Then  els  we  shcdl 
now  show,  any  soliton  solution  will  move  tow'ard  the  back  of  the  pulse,  eventually  falling  into  the 
tail  region  of  A^(t)  and  vanishing.  To  show  this,  just  look  at  the  motion  of  the  center  of  the  soliton, 
which  is  at  Z  =  0.  From  (2)  and  (5),  upon  setting  Z  =  0  and  diff^erentiating,  one  obtains 


dz 


A^ 


dt  ^A^  +  ifu 


<  u 


(7) 


Note  that  this  velocity  is  always  less  than  u,  the  group  velocity  of  the  electromagnetic  waves.  Thus 
the  center  of  the  soliton  always  moves  slower  than  the  wave  envelopes  which  move  at  the  group 
velocity,  u.  Consequently  the  soliton  will  always  move  to  the  back  of  the  pulse  and  must  eventually 
fall  into  the  taul  region  of  the  pulse.  In  the  tail  region,  A{t)  becomes  smaller,  and  the  soltion 
therefore  also  moves  even  slower.  If  A{t)  vanishes,  then  the  soliton  becomes  stationary  zmd  by  (4), 
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the  envelope  amplitudes  of  the  sohton  become  zero.  However,  note  that  even  if  the  field  amphtudes 
vanish,  one  still  will  have  the  medium  excitation  [see  Eq.  (4c)]  nonzero.  What  has  happened  is  that 
the  energy  of  the  soliton  has  been  deposited  in  the  medium  by  leaving  some  of  the  molecules  excited 
(X  ^  0).  And  the  position  of  these  excited  molecules  will  become  stationary.  As  long  as  A(t)  is 
nonzero,  by  (5),  we  have  that  the  SRS  soliton  (and  the  material  excitations)  will  travel  forward, 
albeit  slower  than  the  pump  and  Stokes.  However  as  one  approaches  the  tail  region  of  A{t)  where 
A{t)  is  approaching  zero,  the  soliton  slows  down  [see  Eq.  (7)],  eventually  becoming  stationary  in  the 
laboratory  frame  when  t  is  such  that  A{t)  =  0.  At  the  saune  time  the  electromagnetic  components 
of  the  soliton,  Ai  and  Aj,  have  vanished. 

Of  course  this  brings  to  mind  another  way  to  create  SRS  sohtons.  Namely,  create  an  excitation 
in  the  medium  where  the  area  of  X  is  at  least  7r/2.  Now  before  this  excitation  can  decay,  inject 
the  pump  and/or  a  Stokes  pulse.  One  will  then  see  this  soliton  emerging  on  the  backgroimd  of  the 
Stokes.  This  is  also  easily  predicted  from  the  analogy  between  SRS  and  SIT,  presented  by  Dr.  H. 
Steudel^^  in  these  same  proceedings. 

Let  us  now  consider  some  more  of  the  consequences  of  (1).  As  we  have  already  remarked,  the 
total  photon  number 


A^  —  AJAi  +  A2A2  (8) 

is  conserved,  not  only  for  the  soliton,  but  also  for  the  general  solution.  ^From  (la),  (lb)  and  (8),  it 
follows  that 

d^A  =  0  A  =  A(t)  (9) 

and  A  will  be  only  a  function  of  r.  Then  interpretating  A^  as  the  total  photon  number  gives  us  the 
conservation  of  photon  number  as  the  pulses  propagate  through  a  Raman  active  medium. 

Another  point  is  that  the  evolution  rate  scales  as  the  total  photon  number.  Thus  is  we  double 
the  number  of  photons,  the  interaction  proceeds  twice  as  fast.  This  follows  directly  from  (Ic).  In 
general  we  may  scale  Ai  and  A2  by  A(r).  In  Eqs.  (la)  and  (lb),  since  d^^A  =  0,  these  factors  of 
A(t)  just  factor  out,  leaving  the  equations  form  invariant.  However  in  (Ic)  they  give  (for  7  =  0) 

a,X  =  A^r)l(A,/A){A,/Ar]  (10) 

Now-  define  a  new-  (nonlinear)  time,  T  by 
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(11) 


T  =  /  A^{T)dT 

J-^OC 

then  (10)  becomes 


drX  =  (AifA)(A2fAy  (12) 

which  is  now  independent  of  the  total  photon  number.  So  the  only  effect  of  A{t)  is  to  scale  the 
evolution  rate  in  (10).  It  has  no  effect  on  (la)  or  (lb).  And  if  we  use  the  nonlinear  time  T  to 
’’clock”  the  evolution  as  in  (12),  then  aU  effects  except  two  become  independent  of  A{t).  (Let’s 
also  note  that  back  in  Eq.  (5),  the  last  term  in  Z  is  exactly  this  nonlinear  time,  T.)  Now  these  two 
effects  that  zire  dependent  on  A{t)  are  i)  the  "clock”  rate  as  discussed  above  and  ii)  how  "long”  the 
interaction  will  continue.  To  see  this  latter,  first  note  that  there  is  an  upper  limit  to  T  given  by 


Toe  =  r  A'^{s)ds  (13) 

J^OC 

Now  let  us  only  consider  the  case  where  the  pulse  profiles  are  finite  in  extent  and  contain  a  finite 
number  of  photons.  Consider  the  soliton  solution,  (4),  when  this  Tgc  is  finite.  Then  there  will  be 
a  finite  value  of  x(=  z),  beyond  which  the  soliton  wiU  not  be  foimd.  This  value  is  dependent  on 
the  total  number  of  photons.  If  we  double  the  number  of  photons,  then  T^c  is  doubled  and  the 
soliton  will  double  the  distance  that  it  will  propagate  before  stopping.  The  value  of  Toe  (along 
with  the  amplitude  tj)  determines  the  "clock”  limit  of  the  interaction.  As  soon  as  the  nonlinear 
time  T  reaches  this  limit  of  Toe,  the  interaction  simply  ceases  even  though  the  laboratory  time  is 
continuing.  This  is  true  of  the  general  solution  as  well.  Let  us  represent  the  pulse  components  by 


Then  (1)  becomes 


Ai  =  A(r)e’“cos(i^) 

A2  =  A(r)e’“e“’^sin(^;3) 

a;^[e’“cos(^/3)]  =  -A’c’“c"‘®sin(i^) 
sin(|/?)]  =  X*e’°  cos(ii3) 


(14a) 

(146) 


(15a) 

(156) 


6 


=  ^e^siD^  (15c) 

Now  the  equations  are  clearly  independent  of  A{t)  except  through  the  nonhnear  time  T.  The 
initial  value  problem  to  be  solved  in  SRS  is  given  q,/3  and  0  at  x  =  0  and  X  for  r  — >  — oo,  construct 
the  solution  for  Ai,  A2  and  X.  Integration  of  (15)  will  give  this  solution.  Note  that  this  solution 
will  be  the  same  for  any  profile  A{t)  as  long  as  these  initial  values  are  the  same.  The  only  manner 
in  which  A{t)  affects  the  solution  is  i)  how  fast  the  "ticks”  of  the  nonlinear  clock  are  compared  to 
"ticks”  of  the  laboratory  clock  [AT  =  A'^{T)^t\  and  ii)  the  maximum  time  that  the  nonlinear  clock 
is  to  run  (T  <  T^)- 

The  decomposition  in  (14)  also  illustrates  that  the  d3Tiamical  part  of  SRS  is  simply  the  dynzimics 
of  the  phases;  q,/3  and  6.  The  overall  amphtude  A{t)  scales  out  and  simply  serves  as  a  background 
on  which  the  pheises  evolve.  Of  course,  it  is  necessary  to  measure  to  determine  the  nonlinear 
clock  rate  and  hmit.  But  outside  of  this,  there  is  no  other  need  of  Of  these  phases,  the  most 
important  one  appears  to  be  (3  which  determines  the  relative  ratio  of  the  Stokes  to  the  pump.  In 
most  experiments,  initially  (3  is  very  small  and  as  the  pump  converts  into  the  Stokes,  /?  rotates  to 
either  +7r  or  — tt.  The  other  two  phases,  a  and  0,  seem  to  play  a  lesser  role  and  tend  to  be  slaved 
to  the  major  phase,  (3.  However  it  is  also  more  diflBcult  to  measure  these  phases,  since  one  must 
measure  the  actual  phase  of  the  envelopes  in  order  to  obtain  these.  The  most  that  we  can  say  at 
the  moment  about  these  phases  is  that  the  presence  of  gradients  in  these  phases  probably  will  make 
it  more  difficult  for  SRS  solitons  to  form. 

Lastly,  let  me  make  a  generaJ  remark  in  regard  to  the  inverse  scattering  solution  for  the  SRS 
equations,  Eq.  (1)  with  7  =  0  and  try  to  explain  simply  how  the  1ST  works  for  any  integrable 
system.  To  do  this  it  may  be  simplest  to  take  the  original  forms  of  Scott  and  Chu^,  which  are 


+  i(V2  =  Xv2 

(16a) 

V2x  —  tC^2  =  —X*Vx 

(166) 

Vlr  =  ^{A\Ai  -  AIA2)Vx  -  ^AxAIv2 

(17o) 

— i  % 

V2t  =  '^■AIA2Vx  —  —{A\Ax  —  AIA2)V2 

(176) 
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where  ^  is  caUed  the  spectral  parameter.  Note  that  (16)  and  (17)  are  an  overdetermined  set  of 
equations.  One  may  use  (16)  to  construct  the  solutions  for  Vi  and  V2,  given  appropriate  boundary 
values.  Of  course,  one  may  use  (17)  to  do  the  same  thing.  Now  the  key  question  is  can  one  find  a 
common  solution  to  both  (16)  and  (17).  This  now  is  a  question  of  integrability.  The  way  that  it  is 
answered  is  to  cross-differentiation  these  two  sets  of  equations,  requiring  that  d^drV^  —  drd^Vn  =  0 
for  n  =  1  and  2.  In  other  words,  a  common  solution  to  (16)  and  (17)  will  exist  if  and  only  if  this 
cross-differentiation  does  give  zero.  (A  simple  example  of  this  is  that  a  force  field  can  be  given  as 
a  gradient  of  a  potential  if  and  only  if  the  cml  of  the  force  field  vanishes.)  A  simple  exercise  that 
anyone  may  do  is  to  cross-differentiate  (16)  and  (17).  What  you  will  find  is  that  a  common  solution 
^■ill  exist,  for  arbitrary  C>  ff  and  only  if  the  SRS  equations,  (1),  are  true.  This  has  a  very  simple 
and  important  consequence.  If  by  some  meams  we  can  construct  a  common  solution  for  V\  amd  V2 
which  satisfies  both  (16)  and  (17),  for  arbitrairy  C»  then  the  potentials  in  (16)  and  (17)  must  and 
will  satisfy  the  SRS  equations.  Let  us  restate  this  as  follows:  We  solve  the  SRS  equations  by  not 
solving  the  SRS  equations.  Instead,  we  solve  the  linear  equations  (16)  amd  (17).  If  we  find  any 
solution  for  vi  and  t»2  which  satisfies  both  (16)  amd  (17),  then  we  are  guarenteed  that  Ai,A2  and 
X  will  satisfy  (1)  for  7  =  0.  If  there  is  amy  secret  or  mystery  about  the  1ST,  it  is  contained  in  this 
above  statement. 

A  brief  outline  of  how  this  method  of  solution  is  executed  is  as  follows.  One  first  chooses  one 
of  the  above  pair  of  equations  to  be  the  "eigenvaue  problem”.  For  SRS  we  choose  (17)  to  be  this 
eigen\'alue  problem.  Since  this  problem  is  lineau-,  we  can  use  superposition  and  all  the  other  known 
techniques  for  constructing  solutions  to  this  linear  problem.  In  particular,  we  find  it  very  useful  to 
use  scattering  data  to  represent  the  generaJ  solution  of  (17),  for  arbitrary  Ai  and  A2.  Note  that 
there  is  a  ’'spectral  paraimeter”,  This  pairaimeter  serves  ais  a  "label”  or  "index”  of  the  scattering 
data.  Thus  we  express  the  generail  solution  of  (17)  in  terms  of  scattering  data. 

Now  the  last  trick  is  to  determine  how  the  generad  scattering  data  will  evolve  in  x  by  using 
(16).  Once  this  is  done,  then  one  has  a  general  solution  for  vi  amd  V2  which  satisfies  both  lineair 
equations  (16)  and  (17).  Now  one  can  reconstruct  the  potentials  from  the  scattering  data.  We  do 
that.  Since  both  (16)  and  (17)  are  true,  it  therefore  follows  that  the  potentiails  will  satisfy  (1)  amd 
one  has  then  constructed  the  solution  of  the  nonlinear  SRS  equations  (amd  done  so  by  not  solving 
the  nonlinear  equations).  Of  course,  this  has  been  only  a  simple  description  of  how  it  is  done.  There 
au-e  many  technicail  details  in  do  such  amd  I  shall  simply  refer  you  to  the  quoted  literature  for  the 
mathematical  detaiils. 


8 


ACKNOWLEDGEMENTS 

The  author  wishes  to  thank  H.  Steudel,  C.R.  Menyuk,  and  S.  Burtsev  for  valuable  discussions. 
This  research  was  supported  in  part  by  AFOSR  and  ONR. 

References 

[1]  H.  Steudel,  Aim.  Phys.  Mi  188  (1977). 

[2]  F.Y.F.  Chu  and  A.C.  Scott,  Phys.  Rev.  A  12,  2060  (1975). 

[3]  D.J.  Kaup,  Physica  SH,  143-54  (1983). 

[4]  H.  Steudel,  Physica  £D,  155  (1983). 

[5]  D.J.  Kaup,  Physica  12D,  125-134  (1986). 

[6]  C.R.  Menyuk  and  G.  Hilfer,  Opt.  Lett.  14,  227  (1986). 

[7]  C.R.  Menjoik,  Phys.  Rev.  Lett.  g2,  2937  (1989). 

[8]  G.  Hilfer  and  C.R.  Menyuk,  J.  Opt.  Soc.  Am.  B7,  738  (1990). 

[9]  K.  Druhl,  R.G.  Wenzel  and  J.L.  Carlsten,  Phys.  Rev.  Lett.  51,  1171  (1983). 

[10]  D.J.  Kaup,  A.  Rieman  and  A.  Bers,  Rev.  Mod.  Phys.  51,  275  (1979). 

[11]  H.  Steudel,  (eirticle  in  these  same  proceedings). 


9 


